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A number of ideas and questions related to the construction of quantum pro- 
cesses are discussed. Quantum state extension, entanglement and asymptotic 
(Zh ' behaviour of the entropy are some of the issues explored. These topics are 

C^ ' studied in more detail for a class of quantum processes known as finitely cor- 

^ I related states. Several examples of such processes are presented, specifically a 

'^ '• FVee Fermionic model. 
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[■*-« , 1. Introduction 

rfs ■ The aim of this note is to present a number of ideas and questions related 

.^ , to the construction of quantum processes. For technical convenience we 

f^ ' restrict ourselves to systems with a finite configuration space in a classical 

^3 , context or with d accessible levels in the quantum. 

A reversible dynamics of such an isolated system is rather boring: a finite 
classical configuration space only supports jumps in discrete time while the 
evolution of a finite level quantum systems is almost periodic in time 

id 

Here, the summation runs over the Bohr frequencies of the system, i.e., the 
spacings between the energy levels of the Hamiltonian and the A{uj) are 
the Fourier coefficients of the observable A. 

The evolution of systems weakly perturbed by an environment can be 
reasonably described by a stochastic dynamics, even by a stochastic map F 
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if we only observe the system at regular time intervals. Such maps will typ- 
ically turn pure states into mixed states, a clear signature of their random- 
izing character. In the classical context we deal with matrices of transition 
probabilities and in the quantum setting with completely positive trace pre- 
serving maps. Repeated action of the map on the states yields a dynamics 
{r" I n G N}, Markovian in time. A generic F has a unique invariant state 
Po into which any initial state is driven 

r" (initial state) — > p^. 

Several figures of merit can be defined for this kind of evolution. Two 
well-known ones are the spectral gap 7 of F that controls the asymptotic 
rate of convergence towards the invariant state 

l|r"(p)- Poll ^(1-7)", "large 
and the minimal output entropy 

H'"'"(F) :=min({H(r(p)) | p state} 

Here H denotes either the Shannon entropy of a probability vector or the 
von Neumann entropy of a density matrix. 

A stochastic process in discrete time is a different object: it specifies 
joint probabilities at different times. More precisely such a process is a 
state u; on a half-chain N. At each site of this chain sits a copy of our 
classical configuration space or of the quantum d-level system. We assume, 
moreover, that the chain is stationary: the state u; is invariant under a right 
shift. 

In the classical case the process models a source that is emitting 
every time unit a letter belonging to the alphabet {ei,e2, ■ ■ ■ ,ed}- The 
state a; specifies the probability that the source emits a given word. Let 
(e(0), e(l), e(2), . . .) be a random letter string emitted by the source at times 
i = 0,1,2,... then 

Probf e(0) = eq & e(l) = ei & ... & e(n) = e„j = a;(eo, ei, . . . , e„). 

In the quantum context, the restriction of a; to the first n sites of the 
chain is a density matrix on (8'"A^d- This density matrix encodes all the 
statistical information that can be obtained by applying repeated measure- 
ments on a sequence of n particles emitted by the source. 

The amount of randomness in the process up to time n is quantified by 
the entropy 



For stationary processes H„ satisfies sub-additivity Hm+n < H„j + H„. This 
implies the existence of the average entropy 

h := lim ;^ H„. (1) 

The quantity h is also known as the dynamical entropy of the shift on the 
half-chain, it takes values in [0, logd]. It is a relevant measure of the random- 
ness of the source as, under suitable ergodicity conditions on a;, length-n 
messages can be reliably encoded in a subspace of dimension exp(nh) in- 
stead of d". 

Using strong sub-additivity, H^+m+n + Hm < H^+m + Hm+n, more can be 
proven: the local entropy n n- H„ is increasing, while the entropy increment 
n I— >■ (H„ — H„_i) is decreasing. Both properties fail for general non shift- 
invariant quantum states. As a consequence 

h= lim (H„-H„_i). (2) 

This means that h is not only the compression rate of long messages but also 
the asymptotic entropy production of the source. The importance of (2) is 
that it can be used as a starting point for computing h. This happens for 
some Markov-like constructions where a simple transfer matrix-like con- 
struction generates the n-steps marginal of u) from the (n — l)-th. E.g., 
Blackwell^ has described a procedure for computing the entropy of hidden 
Markov processes and we shall show that an analogous procedure applies 
to free Fermionic processes. 

There are several routes that lead to classical Markov processes, like 
extending two-party states or generating the process in terms of a stochastic 
matrix. We remind here a number of problems and results that arise in 
this context with quantum processes. We also present a general scheme 
for generating quantum processes in terms of a quantum operation on a d- 
level system. Several examples are considered, in particular a free Fermionic 
version. 



2. Classical Markov Processes 

The configuration space of a classical register with d states is just a finite 
set fi — {1,2, ... ,d}. The states are length-d probability vectors 

^l^{^lil),f^{2),...,fi{d)}, fiie)>0, ^^(e) = l. 



The state spaee is a simplex and the extreme points are the Dirac measures 
that assign a probabihty 1 to a configuration. The Shannon entropy 

H(/x):=-5]A^(6)logMe) 

e 

quantifies the randomness in the state. It is easily seen that H is a concave 
function on the state space: 

H(Ai/ii+A2/x2) > AiH(/Xi) + A2H(/x2), A, > 0, Ai + A2 = 1. 

Restricting a state fii2 on ^ composite system 51i2 — fli x 0,2 to the 
subsystem f^i returns the first marginal of fj,i2 



Mi(«i) := 51/^12(^1' ^2) 



The Shannon entropy behaves well with respect to restrictions: 

• monotonicity: ^{fJ^i) < H(/X]^2): 

• sub-additivity: H(/X]^2) ^ ^ifJ-i) + H(/X2), and 

• strong sub-additivity: H(/Xj23) + ^(^2) ^ ^(/Xjj) + ^{^''23)- 

We can now consider the following state extension problem. Suppose 
that we are given two probability vectors fj,i2 ^-iid i'23 that agree on the 
middle system: /Xj = 1^2- Can we find a joint extension for /Xj^j ^^'^ 1^23? 
More explicitly: can we find a state ^]^23 ^^ ^123 that restricts to (1^2 
on rii2 and to 1/23 on ri23? This is indeed possible and clearly the set of 
joint extensions ^^^23 i^ convex and compact. We can therefore refine the 
question and ask for a joint extension of maximal entropy. A straightforward 
computation yields the answer: 

/ N Mi2(ei>e2)i^23(e2,e3) Mi2(ei, ^2) «^23(e2, £3) ,^. 

Pi23(ei> £2, £3) :- — 7-T — 7-T • (3) 

fJ'2\^2) J^2(e2J 

Actually, P123 saturates the strong sub-additivity inequality: 

H(P123) + H(P2) = H(P12) + H(P23)- 

Unsurprisingly, there is a direct connection with thermal equilibrium states. 
If we introduce Hamiltonians 

P12 = c"''l^ 1^23 = e"'''% and /X2 = 1^2 = c"''% 
then 

„ _ „-('il2+'i23-'l2) 

P123 ~ 6 



Let us start with a two-party probability vector /x that is shift-invariant: 

We can repeatedly apply the Markov extension procedure (3) to get a sta- 
tionary process 

a;(eo,ei,...,e„) = . ^ . ^ 7 ^ • (5) 

Another procedure is to start with a d x d stochastic matrix T. The 
entry T^jg^ is the probability for jumping from state ei to 62, therefore 

T,,,, >0 and ^Te.e^ = 1- (6) 

The invariant state /x is a row vector determined hy fiT — fi. The Markov 
process is now obtained by putting 

a;(eo, ei, •••,£«) =M(eo)r'eoei •••T'e„_i£„- (7) 

Both constructions (5) and (7) agree if we put 

rp _ M(ei,e2) 
/x(ei) 

The rows of a stochastic matrix T are probability vectors. The minimal 
output entropy of T is simply 

H'"'"(r) = smallest entropy of rows of T 

while the entropy of the process is a smooth version of this quantity 

h — /x-average of entropies of rows of T. 

3. Extending Quantum States 

When turning to quantum state extension the situation gets more compli- 
cated. Quantum states allow for more freedom, as they exhibit correlations 
that are not present in classical systems, but this imposes at the same time 
more stringent positivity conditions. 

States on a full matrix algebra Aid can be identified with density ma- 
trices: non-negative matrices with trace one. The convex set of density 
matrices is very unlike a simplex. A density matrix that is not an extreme 
point of the state space, i.e., that is not a one-dimensional projector, allows 
many decompositions in extreme states. In contrast with classical systems 
such a state can therefore not be seen as a well-defined ensemble op pure 



states. We need d^ — 1 real parameters to describe the state space of Md 
while 2d— 2 parameters suffice to label the pure states. This means that the 
boundary of the state space contains many flat subsets. Nevertheless the 
pure states form a very nice smooth manifold. The case of a single qubit is 
exceptional: its state space is affinely isomorphic to the Bloch ball by the 
standard parametrization 

p = i (1 + X • cr), X e ]R^ ||x|l < 1. (8) 

In this case, every point of the boundary is also an extreme point. For higher 
d, a smooth parametrization of the pure states does not define a boundary 
of a convex set. 

For a composite system, restricting to a sub-system amounts to taking 
partial traces over remaining parties 

pi := Tr2/9i2. 

The entropy of a state with density matrix p is given by the von Neumann 
entropy 

H(p) = -Trplogp. 

However, already the most basic property of the Shannon entropy, mono- 
tonicity, does not carry over. Consider for example the maximally entan- 
gled two-qubit state |$+)($+| with |$+) = (|00) + |ll))/y2. Its entropy 
H(|$"'")(<I>"'"|) is zero as it is a pure state, while its restriction pi is the max- 
imally mixed state, which has maximal entropy, so clearly H(pi) ^ H{pi2). 

An important property that holds both for classical and quantum sys- 
tems is that if the marginal pi of a bipartite state is pure then pi2 = pi®p2- 
This is an important ingredient of the theory: it namely allows to isolate a 
system from the rest of the universe. At the same time it is also a severe 
constraint on quantum systems because there are plenty of pure states of a 
composite system. In particular the restriction of an entangled pure state 
can never be pure and we can therefore not separate a party of an entangled 
system from the outside world, which is more or less what goes wrong with 
the locality assumption in the EPR paradox. 

Factorisation of extensions of pure states has also a bearing on joint 
extensions of states as considered in the previous section.^ Indeed, suppose 
that pi2 and p23 are pure and agree on the middle system, which is easily 
feasible, then a joint extension /3123 can only exist for pi2 and P23 pure 
product states. Therefore a generic pure two-party state with inner shift- 
invariance as in (4) cannot be extended. 



Suppose that density matrices pi2 and 1723 agree on the middle system 
and can be jointly extended. The set of extensions is still convex and com- 
pact and so we can still look for the maximal entropy extension. Finding 
this state is hard, however, because generally 

[pi2 (» I3 , ll (» cr23] 7^ 

or, equivalently, if pi2 and (T23 are equilibrium states corresponding to 
Hamiltonians /112 and /123 

Tr3 exp(/ii2 + /123) 7^ cxp(/ii2 + /i2) . 

Moreover, the maximal entropy extension will not saturate the strong sub- 
additivity. 

Actually, a nice characterisation of equality in strong sub-additivity for a 
state /9i23 on a space 'Hi®'H2®'-H3 in terms of decompositions of the middle 
space has been obtained in the paper by Hayden et al.^. The necessary and 
sufficient condition is that the middle Hilbert space 'H2 decomposes as 

H2 - e Hfeft <» Height and P123 - e A„ p?2 ^ P23 

Q a 

with {Aq} convex weights. 

A Qubit Example with SU (2) -symmetry 

An example of the limitations imposed on quantum state extensions can be 
worked out for qubits with a SU(2)-symmetry. In order to impose SU(2)- 
symmetry on single qubit observables we use the adjoint representation of 

SU(2) 

Ad{U):A^^UAU*, UeS\J{2), A e M2. 

This is a reducible representation that decomposes into a spin and a spin 

1 irrep: 

M2 = Cl®Ccr. 
The only SU(2)-invariant state on AI2 is the uniform state 

For 2 qubits Ad{U (E) U) decomposes into 2 spin 0, 3 spin 1 and 1 spin 

2 irrep. There exists now a one-parameter family of SU(2)-invariant states 

p=i(l-A)(l-p) + Ap, 0<A<1. 



Here p is the projector on the singlet vector -y= (|10) — |01)) in (D^ (8) (D^. 
This projector commutes with every unitary of the form U ®U and every 
two-qubit observable that is SU(2)-invariant is a linear combination of p 
and 1. Clearly, SU(2)-invariant two-qubit states satisfy 

< (p> = A < 1. 

The two-qubit state p is separable for < A < ^ and entangled for i < 
A < 1. Hence the expectation value of this projector for a certain process 
tells us how much bipartite entanglement between two neighbouring spins 
is attainable. 

For 3 qubits the SU(2)-invariant states can still easily be determined 
but things become more complicated with increasing number of parties. 
Let 

pi = p (g) 1 and p2 = 1 (g) p and put 

q = I (Pl + P2 - PlP2 - P2Pl)- 

The algebra of three-qubit obscrvables that are SU(2)-invariant is not 
Abelian. It can be decomposed into a direct sum of (D and M.2 where 
(D is identified with (Dq and M2 with the algebra generated by pi and p2, 
not including 1. An SU(2)-invariant three-qubit state is of the form 

p=\{l- A)(l -q) + A(api + 6p2 + cpip2 + cp2pi) 

with 

0<A<1, a,6eK., ceC, 2a + 26 + 3fJe(c) = 1, tmd \c\^ < Aab. 

If wc look for a SU(2)-invariant three-party state with partial shift- 
invariance, then we find the following constraint on the expectation of p 

< (pi> = (P2> < 3/4. (9) 

SU(2) and shift-invariant states on more parties will satisfy stronger 
upper bounds on the expectations of p, see (9). Ultimately, if we look 
for a shift-invariant extension on the full half-chain then, using the Bethe 
Ansatz*, one can show that 

0< (p) <log2«0.69. (10) 

We may look for the largest expectation value of p that can be obtained 
within classes of shift-invariant states that can easily be handled. Consider 
as a first example point-wise limits of shift-invariant product states. Such 
states are actually invariant under arbitrary finite permutations of sites 



on the half-chain and are usually called exchangeable. Using the Bloch 
parametrization (8) and 

p = i(l — CTi • (T2) with <Ti = <T (g) 1 and (T2 = 1 (g) cr 

we have to maximize 

xeR^K^ ^ Tr[(l+x-(Ti)(l+x-cr2)(l-cri-cr2)] 

subject to the constraint ||x|| < 1. It is easily seen that the maximum is 
reached for x = for which value (p) = ^ Hence 

(p) < J for exchangeable states. (11) 

The largest expectation for p that can be reached within the class of 
product states is 

max(^^ Tr[(l + xi ■ (Ti){1 + X2 ■ (T2)(1 - cti -era)] 

subject to the constraint ||xi||, ||x2|j < 1. The maximum i is attained for 
Xi = — X2 = X where X € IR'^ is an arbitrary vector of length 1. Therefore 

(p) < ^ for separable states. 

Moreover, this maximum is attained for shift-invariant separable states that 
are equal weight mixtures of period-2 product states 

i |eo)(eo|®|ei)(ei|«)|eo)(eo|®- • -4-^ |ei)(ei|® |eo)(eo|® |ei)(ei|®- • • , (12) 

where {eo, ei} is any orthonormal basis in (D^. Hence 

(p) < 2 ^'^^ shift-invariant separable states 

is an optimal upper bound. States of the form (12) are extreme shift- 
invariant states which allow a convex decomposition in clustering period-2 
states. This is called Neel order of period 2. The value ^ for shift-invariant 
separable states is still not close to the maximum value of log 2. One can 
get closer by constructing more general quantum processes. 

4. Constructing Processes 

We now turn to the construction of classical and quantum processes using as 
initial data a unity preserving CP map T : Aid ^^ Md with invariant state 
p. In the classical case this reduces to a stochastic matrix T with invariant 
measure fi. The construction is based on finitely correlated states^, also 
called matrix product states. These are processes for which the correlations 
across any link can be modelled by a finite dimensional vector space. These 



10 

states are more general than the ones we have considered until now and are 
easily constructible in the thermodynamic limit, unlike the Bethe Ansatz 
states. The finitely correlated states where e.g. used in the lectures by 
J. Eisert under the form of tensor networks. Actually pure states have 
been considered there as these lectures were focusing on ground states. We 
present here a different version that is adapted to mixed states. 

The starting point is a unity preserving completely positive (UPCP) 
map 

that is compatible with the given F in the following sense 

A{A(E)1)^ A{l(g> A) = Ti A), A eMd- (13) 

A process a; is now generated by repeatedly contracting the local observ- 
ables on the half-chain. Consider a sequence of UPCP maps 

A(o) ■.= A:Md®Md^Md 

A^^' :== A o (A «) id) : Md «> {Md «> Md) -^ Md 

■ (14) 

A(") :^ a o (A("-i) ® id) : Md ® { Md®- - -®Md) ^ Md- 

(n + 1) times 

The expectation of a local observable A„ € ®^M.d is then computed as 

u;(A„) := Tr{p A(")(l ® A„)}. (15) 

To define a stationary process, (15) must satisfy a number of re- 
quirements. The definition should be consistent in the first place, namely 
u){An (E) 1) = u){An)- This follows from the compatibility (13) and the 
invariance of p: 

uj(A„ (E)l)= Tr|p A("+i) (1 ® A„ ® 1)} 

= Tr{p (a o (A(") ® id)) (1 «) ^„ ® 1)} 

= Tr{pA(A(")(l®^„)®l)} 

= Tr{pr(A(")(l®yl„))} 

= Tr{pA(")(l®AO} 
= w(A„). 



11 

Next, we need positivity. This follows immediately from the complete pos- 
itivity of A. The com.patibility condition im.plies that A m.aps the identity 
on (C^ (E) (D"^ to the identity on (D''. This implies the normalization and 
stationarity of u). 

It is important to observe that compatibility (13) imposes a severe re- 
striction on r. Not every UPCP transformation F of Aid admits a com- 
patible extension. Moreover, the compatible extensions of F, whenever such 
extensions exist, form a compact and convex set and one may wonder about 
particular extensions. We now turn to some classes of examples. 

4.1. Hidden Markov Processes 

A classical observable, i.e., a R- valued function / on configuration 
space n — {1,2, ... ,d} is naturally tabulated into a vector f ~ 
(/(I), /(2), . . . , /(d)) e M!^ and identified with a diagonal matrix in Aid 
through the map 



dia(f)=5^/(6)|e)(6| 



The relation between a (completely) positive transformation F of Aid and 
a stochastic d x d matrix is then 

F(dia(f)) =dia(rf). 

This allows to rewrite the compatibility equation (13): a stochastic matrix 
S* : R'' «) R'' ^ R'' is compatible with a stochastic matrix T : R'' -> R'' if 

Let us introduce d square matrices of dimension d with non-negative entries 

E{e)ip^ri — '-'(^,(jj,£)- 

The process generated by S is then seen to be 

a;(eo, ei, . . . , e„) = (/x , E{eo)E{ei) ■ ■ ■ E{e„)l), 

where 1 € R'' has all its entries equal to one and //, is the invariant proba- 
bility vector for T. 

A stochastic matrix T always allows the extension 
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The corresponding process is the usual Markov process (7). More general 
extensions a; are hidden Markov processes: there exists a larger configura- 
tion space fti, a function F : ili ^S' H, and a Markov process Wi on fti such 
that 

w(eo,ei,...,e„) = ^ uJi{(po,ifi,...,ifin)- 

The entropy of hidden Markov processes can be computed using a method 
due to Blackwell^'^. The starting point is the asymptotic entropy production 
formula (2). The construction of the process, adding one point at a time, 
see (14) and (15), defines a dynamical system on the length-d probability 
vectors. The entropy of the process is then obtained as an average over 
entropies of probability vectors with respect to the invariant measure of the 
dynamical system. Numerical evidence suggests that the Markov extension 
has the smallest entropy amongst all. 

4.2. Qubits with SU(2)-invariance cont. 

In order to have manifest SU(2)-invariance of the process we impose SU(2)- 
covariance both on the CP transformation of A^2 and on its compatible 
extensions from Al2<8'A^2 to 7^2- Let Q ^^ Ughe a. unitary representation 
of a group Q on s. Hilbert space %. The adjoint representation lifts it to a 
representation of Q on the bounded linear transformations B{'H) of H: 

AdiUg){A) - UgAu;, geg, Ae B{n). 

Given two unitary representations U^^' and C/"' of Q on "Hi and 'H2 a map 
r : B{Hi) -^ B{H2) is covariant if 

roAd(t/(i)) = Ad(c/(2))or (16) 

The Choi-Jamiolkowski encoding of a linear map F : Alc(j ^ Al^j is 
very convenient for handling complete positivity 



C(^):=^|^)(J|®^(|^)( 



r is completely positive if and only if C(r) is positive semi-definite. The 
encoding depends on the chosen basis through the matrix units \i){j\ but 
only up to unitary equivalence as 

C(r o Ad(C/)) = Ad(t/"^ ® 1) o C(r) and 

c(Ad(t/) o r) = Ad(i (E)U)o c(r). 
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The covariancc condition (16) for F : Aidi -^ -Md-j translates for its Choi- 
Jamiolkowski encoding into 

[ul'^ ® c/f ) , c(r)] ^0, geg. 

Here A is the complex conjugate of the matrix A. For SU(2) there is an ad- 
ditional simplification because the conjugate of SU(2) is unitarily equivalent 
to SU(2). 

It turns out that there is a one-parameter family of SU(2)-covariant 
UPCP transformations of AI2 

F(cr) =^cr, -i < Ai < 1- 

The SU(2)-covariant UPCP maps A : M2 ® M2 -^ M2 compatible with F 
are parametrized by three real parameters 

A(<Ti • (T2) ^ at, 

A(cri) = A(cr2) =F(cr) = ^<T, and (17) 

A(<Ti X (T2) — rjcr. 
Complete positivity imposes constraints on a, /i, and 77 

le^i - a| < 3 and 3 - 2a - a^ + 12^ - 12a^ - Qrf > 0. (18) 

These conditions can be obtained by imposing positivity on the Choi matrix 
of A. The allowed region is a piece of a cone in IR'^. We then compute the 
expectation of p 

(P> = i - i (-^i • ^^2) - i - i Tr ^ A(cT-' A(cT-' 1)) 

7 

The maximum in the allowed parameter region is attained for a = — | and 
IJ' — 1 and is independent of 77. Therefore 

(P) < i (19) 

for ( ) a stationary and SU(2)-invariant process as in (14). This should be 
compared with (11). 

In passing from exchangeable to shift-invariant separable states we actu- 
ally allowed product states of period 2. This can also be applied to processes 
of the type (17). Considering <ti • <T2 as the contribution to the energy of 
two neighbouring spins, a minimal value of {cri ■ CT2) corresponds to a maxi- 
mal value of (p) and this is expected to happen for spins as anti-parallel as 
possible. Therefore the second requirement in (17) is inappropriate and we 
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should consider general SU(2)-covariant maps A : A^2 ® M.2 -^ M.2- These 
are determined by four real parameters 

A(cri • (T2) = al, 

A(cri)=^cr, h{cT2) = vcr, and (20) 

A(cri X (T2) = rjCT. 

Complete positivity imposes the constraints 

|3u + 3i^ — Qfl < 3 and 

(21) 
3 - 2q; - a^ + 6(1 - a)(Ai + z^) - 9(m - i^)^ - 977^ > 0. 

Wc now introduce two SU(2)-covariant maps A^ : M.2 ® M.2 -^ M.2 as 
in (20) with defining parameters (a^, /i^, Vi, rji). The expectation of p in the 
equal weight average of these period-2 processes is given by 

(P> = i - i ("2M1 +aii^2)- 

Maximizing this in the allowed parameter region yields 

(p) = I = 0.625 

which is within 10% of the optimal bound and well within the entangled 
shift-invariant states. 

4.3. Davies Maps 

An interesting and physically relevant class of channels are the Davies maps, 
they arise in the reduced dynamical description of a system with a discrete 
level structure weakly coupled to a thermal bath^. The level structure of the 
small system is preserved in the sense that such a map F is parametrized by 
a stochastic map T and a decoherence matrix D. The matrix T describes the 
stochastic evolution of the diagonal elements while D gives the damping of 
the off-diagonal terms. Assuming that the system Hamiltonian is diagonal 
in the canonical basis 

r(dia(93)) =dia(r93) and r(ey) = Aj-ejj, i^j. (22) 

Here, e^ = |ei)(ej|. Moreover, T is detailed balance and D is real symmet- 
ric. Detailed balance means that T is Hermitian for the stationary measure 
fjb that is interpreted as the Gibbs state of the system 

tiifTig)) = /x(r(/).9), / and g real-valued. 

This condition is equivalent with micro-reversibility: the occupation rate of 
level i times the jump probability from i to j is equal to the occupation rate 
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of j times the jump probability from j to i. Another equivalent condition 
is 

TijTjkTkt = TikTkjTji for aU choices of i,j, k. 

Complete positivity additionally imposes that 



D2I T22 ■ ■ ■ D2d 



\DdiDd2 ■■■ TddJ 



be positive semi-definite. (23) 



The action of F is quite clear: decay of the off-diagonal elements and birth 
and death process for the diagonals. The relations between the different 
rates are encoded in the positivity condition (23). E.g., one can readily 
check that the decay rate of the off-diagonals cannot be less than half the 
rate of convergence to equilibrium for the diagonal process. 

For Davies maps one could expect the standard basis vectors to be 
minimizers of output entropy but this is not generally true. It has been 
shown* that already for a single qubit a true superposition of ground and 
excited state is the minimizer in a regime where the map is close to the 
identity map and so truly quantum. High powers of a Davies map converge 
to the projector on the equilibrium state which is entanglement breaking. 
In this regime the minimizer for output entropy is the state corresponding 
to the row of minimal entropy in T. 

The construction of a process as in (14) requires a Davies map rather 
closer to the projector on the equilibrium state than to the identity. For a 
single qubit 

T = P , " ° j , with < 6 < a < 1 

and with d the damping factor of the off-diagonal element one checks that 
it generates a process if 

d^ < i(l-a)(l-6). 

4.4. Free Fermionic Processes 

For both Bosons and Fermions there exists a notion of Gaussian states and 
maps^'^*^ that are considerably simpler to handle than general ones. The 
names free, quasi-free, quadratic, linear, and determinantal are also used. 
Moreover, these states and maps are good approximations whenever the 
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statistics dominates over the interactions. A considerable benefit is also the 
scaling behaviour: the dimension of the free objects grows linearly in the 
number of particles instead of exponentially. We shall here only describe 
the defining free objects without connecting them to the true observables, 
states, and maps of a many particle system. This yields a kind of nieta- 
description. Moreover, we restrict ourselves to Fermions. 

The observables of Fermions with mode space H are the trace-class 
operators Ti{'H) on %. Apart from their linear structure commutators are 
also useful. The mode space 'H12 of a bipartite system is just the direct 
sum of the mode spaces "Hi and "^2 of the corresponding subsystems and 
observables of system 1 are extended by putting 

A e IiiUi) H^ A ® e IiiHu). 

The symbols play the role of states, they are operators Q onH such that 
< Q < 1. The expectation of an observable A is just TrQ A. Mixtures of 
symbols are constrained by the following requirement: let < A < 1 then 
the mixture of Qi ^ Q2 with weights A and 1 — A can be formed if and only 
if <5i ^ Q2 is of rank 1 and it yields the symbol AQi + (1 — A)(52- It then 
follows that a symbol is pure if and only if it is an orthogonal projector, 
possibly 0. Given a composite system 'H12 = "Hi © "^2 and two symbols Qx 
and Q2, the product state has symbol Qi (BQ2- In this context, a separable 
state is just block diagonal in the mode space decomposition. We shall also 
need the von Neumann entropy of a symbol on a finite dimensional space 

H(0) :=-Tr[QlogO + (l-Q)log(l-Q)]. (24) 

In particular, H{Q) = if and only if Q is pure. 

A trace-preserving completely positive map from "Hi — > "7^2 is a couple 
{A, B) of linear maps where 

A:Hi^H2, B-.T-Li^Hi, tmd Q < B <1 - A* A. (25) 

The action on a symbol is given by 

Q^Q' = A*QA + B. 

Observe that such maps are compatible with the notion of convex mixture 
of above because there is only a single Kraus-like operator appearing in (25) . 
Composition of free CP maps is given by a semi-direct product 

(A, B) o (A', B') = {AA,B' + [AfB A). 

We now mimic within the Fermionic free context the construction of 
a process starting from a CP transformation [A^B] of (D''. Such a map 
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can be extended to a compatible map (C, D) from C* ® (C^ -^ (C^ if and 
only if A* A < min({i 1, 1 — B}) and the extensions are labelled by an 



C={AA) and D 



B X 

X* B 



The symbol Q invariant under {A, B) is the solution of 

Q^ A*QA + B. 

It is not hard to see that the outcome of the construction is a symbol Qo 
on ©j^C'' which is a block Toeplitz matrix: the dx d matrix entries in Q^ 
are constant along parallels to the main diagonal. Explicitly 



(Qo 



Q and (Qo 



' i i+n 



{A*Y\Q-B + X). 



(26) 



The entropy can now be computed in two different ways, we can either 
compute the limiting average entropy as in (1) or the asymptotic entropy 
production as in (2). The first method relies on an extension of Szego's 
theorem to block Toeplitz matrices. For the second we need either a much 
finer control on the spectra of principal sub-matrices of a block Toeplitz 
matrix, which appears to be hard, or we have to exploit the smoothness of 
the entropy function. We follow this last approach. 

Let T : [— tt, 7r[— > M.d be an £°°-function taking values in the Hermitian 
d X d matrices and put 



T : 



( f(o) f(i) r(2) 
f(-i) f(o) r(i) 
f(-2) r(-i) r(o) 

V ; ; ; 

where T are the Fourier coefficients of T 



\ 



f{k) 



27r 



d0Tie)c 



-2TTikB 



An extension of the classical Szego theorem^ ^ reads 

Theorem 4.1 (Szego). For any continuous complex function f on ]R 



hm - Tr/(P„rP„) = ^ I d9 Trf{T{9)) 



(27) 



OO/f^d 



where P^ projects onto the first n terms in ©^(D 
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The theorem gives us information on the main asymptotic behaviour of 
the eigenvalues of principal sub- matrices of T . Consider e.g., the case d = 1, 
then Szego's theorem can be rewritten as 



lim - Tr /(P„rP„) 



Kdx) f{x) 



M 



where 



/x(] -oo,a;]) = 



27r 



de. 



T(e)<x 



If we order the eigenvalue list {A„.j} of P„TP„, then {A„.j} interlaces 
{K+i.j} and 



■ lim ■ 



>A„,, =M. 



J = l 



A fine asymptotic control on the eigenvalues could be used to obtain the 
average (27) as an asymptotic growth 

lim {Tr /(P„fP„) - Tr/(P„_if P„_i)}. 

Numerical evidence, however, shows that the behaviour of eigenvalue spac- 
ings can become erratic when T oscillates. 

In fig. 1 the function T{e) = \ + \ cos(6') + \ sin(26l) is plotted. In fig. 2 
the eigenvalue lists of the first 50 principal sub-matrices are shown together 
with a plot of the eigenvalues of the 100 x 100 sub-matrix. This last plot 
approximates well the reordered function T. 




Fig. 1. The function T(0) = | -f I cos{e) + \ sin{20) 
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Fig. 2. Eigenvalues of principal Toeplitz sub-matrices 

Strong sub-additivity of the entropy guarantees the existence of the 
asymptotic entropy production but a much more general result can be 
proven by extending Szego's theorem to reasonably smooth functions. This 
result can then be applied to the computation of the entropy of our pro- 
cesses. 

Theorem 4.2. For a block Toeplitz matrix T and an absolutely continuous 
complex function f onTR 

lim (Tr/(F„+iTF„+i)-Tr/(P„fP„)) = -^ / d0TTf{T{e)). 



As a corollary we get 






with 



^joo — ^ 



1 ~ A*e'-' 



d9H{Qoom 



{Q-B + X) + h.c. 



and H as in (24). 



5. Conclusion 

The construction of quantum processes is a lot less straightforward than 
for their classical counterparts. The intricate nature of quantum correla- 
tions complicates even the seemingly simple task of finding extensions of 
overlapping states. 

Processes that can nevertheless be easily constructed, like exchangeable 
or separable states, are not general enough to study interesting quantum 
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behaviour. On the other hand, more general constructions Uke the Bethe 
Ansatz become difficuh to handle as the size of the process increases. 

The processes we have studied here, the finitely correlated states, lie 
somewhere in between the previous two classes. By construction, they are 
well-behaved as the length of the process grows. We have also seen by study- 
ing some concrete examples that such states do in fact exhibit interesting 
quantum characteristics. 
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